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Measurements  on  BaTiO-  ceramic  were  made  of  inverse  dielectric 
constant  versus  uniaxial  strain  for  values  of  strain  ranging  between 
•0,0045  and  -0.015,  The  strain  was  produced  in  the  region  behind  a 
shock  front  generated  In  a  slab  of  the  material  by  impacting  It  with  u 
flying  plate.  The  plate  was  accelerated  in  a  gas  gun.  The  slabs  were 
disks  with  electrodes  on  the  flat  surfaces.  The  normal  stress  produc¬ 
ing  these  strains  ranged  between  8.5  and  26  khars.  The  dielectric  con¬ 
stant  was  obtained  from  oscillograms  of  the  shock -induced  depolariza¬ 
tion  current  throngh  a  circuit  connecting  the  electrodes.  A  straight 
line  was  fitted  to  the  curve  of  Inverse  dielectric  constant  versus 
strain.  The  slope  of  the  line  was -0.11  *  0.04.  Thermodynamic  theory 
was  used  to  calculate  this  slope  assuming  the  material  wus  in  the  para- 
electric  phase  for  the  range  of  experimental  shock  pleasures.  The 
calculated  value  was -0.13.  The  calculated  value  is  in  agreement  with 
the  experimental  value  and  the  assumption  that  the  material  is  in  the 
paraelectric  phase  for  shock  pressures  botweon  8.5  and  26  kbars.  The 
agreement  implies  that  the  transition  at  approximately  room  temperature 
is  produced  by  a  shock  pressure  less  than  or  equal  to  about  8,5  kbars. 
This  is  markedly  different  from  the  experimental  results  for  hydro¬ 
static  pressure  where  the  transition  occurs  at  20  kburs  at  room 
temperature, and  the  material  is  in  its  paraelectric  phase  above  20  kbars. 
In  both  casos  the  temporuture  is  approximately  room  temperature.  The 
transition  strain  and  shock  prossure  was  calculated  using  thermodynamic 
theory  and  assuming  that  the  transition  could  be  approximated  as  a 
socond  order  one.  The  result  was  a  transition  strain  of  -0.0044  and  a 
transition  shock  pressure  of  8.6  kbars.  This  is  consistent  with  the 
preceding  results  and  also  supported  by  experimental  findings  by  Doran 
who  observed  a  discontinuity  in  sound  velocity  at  that  shock  pressure. 
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1. 


INTRODUCTION 


Almost  all  solid  materials  undergo  a  yield  under  sufficiently 
high  shock  compression.  This  yield  manifests  itself  as  a  change  In 
slope  of  the  shock  pressure  versus  compression  curve.  The  pressure 
at  which  yielding  occurs  is  called  the  Hugoniot  Elastic  Limit  (HEL) , 
and  is  characteristic  of  the  particular  material.  When  the  shock 
pressure  is  considerably  above  the  HEL,  the  configuration  of  stress 
and  strain  can  generally  be  assumed  to  approximate  those  produced’  by 
hydrostatic  pressure  (ref  1).  A  polymorphic  solid-solid  transition, 
which  occurs  at  a  shock  pressure  considerably  above  the  HEL,  can 
reasonably  be  expected  also  to  be  produced  by  a  hydrostatic  pressure 
equal  in  magnitude.  This  expectation  has  been  verified  in  several 
instances  (ref  2-4).  Below  the  HEL,  however,  the  strains  behind  the 
shock  front  can  no  longer  be  expected  to  approximate  those  due  to 
hydrostatic  loading.  The  strain  now  is  uniaxial,  involving  no  plastic 
flow  and  no  molecular  rearrangement  (provided  there  is  no  transition). 
The  stresses  are  nov.  determined  by  the  elastic  stiffness  coefficients 
and  the  requirement  that  the  strain  be  uniaxial.  In  general,  this 
means  that  it  would  be  coincidental  if  below  the  HEL  a  phase  transition 
produced  by  a  certain  hydrostatic  pressure  would  also  be  produced  by 
a  shock  pressure  of  the  same  magnitude,  since  the  transition  pressure 
depends  on  the  relative  free  energy  of  the  competing  phases  and  the 
free  energy  depends  on  the  stresses  and  strains. 

The  phase  transition  considered  in  this  report  is  the  ferro- 
electric-paraelectric  transition  in  BaTiO-j.  The  material  is  a  ferro¬ 
electric  with  a  tetragonal  structure  below  approximately  120°C 
(at  atmospheric  pressure)  and  a  paraelectric  material  with  a  cubic 
structure  above  120°C.  The  material  used  in  the  experiments  described 
is  the  polycrystalline  ceramic,  composed  of  relatively  small  crystals 
sintered  together.  The  density  is  within  a  few  percent  of  the  crystal 
density.  Increasing  hydrostatic  pressure  reduces  the  transition 
temperature  of  both  single  crystal  and  ceramic  at  a  rate  of  about  4°C 
per  kbar.  This  was  first  observed  by  Merz  (ref  5)  and  has  been 
thoroughly  investigated  by  Samara  (ref  6).  At  room  temperature,  the 
ferroelectric-paraelectric  transition  occurs  when  a  hydrostatic 
pressure  of  20  kbars  is  applied  to  either  the  ceramic  or  the  single 
crystal  material  (ref  6-10).  The  pressure  at  which  this  transition 
takes  place  and  also  the  slope  of  the  inverse  dielectric  constant  versus 
hydrostatic  pressure  curve  can  be  calculated  from  thermodynamic  theory, 
and  the  results  obtained  are  generally  in  agreement  with:  the  results 
of  the  hydrostatic  measurements. 

In  this  report,  both  experimental  measurements  and  calculations 
are  described  which  are  analogous  to  those  that  have  been  made  for 
hydrostatic  pressure.  This  new  work,  however,  deals  with  the  dependence 
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of  the  f erroeleetric  puraelectric  transition  on  uniaxial  strain  rather 
than  hydrostatic  pressure.  The  experimental  work  described  consists 
of  measurements  of  the  inverse  dielectric  of  the  BaTi03  ceramic  as 
a  function  of  uniaxial  strain.  The  uniaxial  strains  are  produced  by 
shock  waves  created  in  ceramic  disks,  using  the  impact  of  a  flying 
plate  to  produce  the  shock.  To  produce  uniaxial  strain,  the  shock  pres¬ 
sures  are  kept  below  the  HEL  of  the  material,  which  is  approximately 
25  kbars  (ref  11,13).  Temperatures  behind  these  fronts  are  just  slightly 
above  room  temperature  (ref  12). 

The  inverse  dielectric  constant  is  obtained  from  oscillograms  of 
shock-induced  depolarization  currents  through  an  effective  short 
circuit  connecting  electrodes  on  the  ceramic  disk,  Thermodynamic 
calculations  are  made  for  the  slope  of  the  inverse  dielectric  constant 
versus  strain  and  also  of  the  transition  pressure  as  a  function  of 
uniaxial  strain.  The  experimental  results  are  in  approximate  agreement 
with  the  theory.  The  experimental  and  calculated  results  for  the  shock 
pressure  are,  However,  very  different  from  those  for  hydrostatic  pressure, 
and  it  is  quite  obvious  that  for  the  case  of  uniaxial  strain,  shock 
pressure  and  hydrostatic  pressure,  in  general,  are  not  equivalent  as 
far  as  transitions  are  concerned. 

There  is  a  practical  design  reason  for  knowing  at  what  shock 
pressure  f erroelectric-paraelectric  transition  occurs.  Barium  titanate 
ceramics  are  used  as  piezoelectric  power  sources  for  contact  fuzing. 

The  power  is  developed  as  a  depolarization  current  produced  by  an 
impact-created  shock  wave  in  ceramic.  The  f erroelectric-paraelectric 
transition  produces  complete  depolarization  in  the  ceramic.  Knowledge 
of  the  shock  pressure  at  which  the  transition  occurs  is  an  important 
asset  when  designing  point-contact  fuzes,  since  the  shock  pressure 
depends  on  the  impact  velocity  and  the  impacting  materials.  Knowl¬ 
edge  of  this  pressure,  which  produces  full  depolarization  of  the 
ceramic,  is  even  more  important  than  might  first  appear.  There  are 
experiments  which  indicate  that  for  strong  shocks  there  is  marked 
increase  in  the  conductivity  of  the  ceramic  element,  shorting  it,  and 
thus  reducing  the  charge  released  by  depolarization  into  an  external 
load  (ref  14).  The  shock-transition  pressure,  therefore,  is  not  only 
the  minimum-pressure  for  maximum-charge  release  but  also  a  pressure 
above  which  charge  release  may  be  expected  to  decrease.  As  a  result 
maximum  charge  release  is  produced  over  a  limited  pressure  rant?*. 


The  next  section  describes  and  develops  the  thermodynamic  theory 
relevant  to  this  report. 


2. 


THERMODYNAMIC  THEORY 


In  a  thermodynamic  theory,  the  free  energy  is  expressed  as  a 
function  of  independent  thermodynamic  variables.  Devonshire  (ref  15) 
has  expressed  the  Helmholtz's  free  energy  for  the  case  of  cubic 
symmetry  in  terms  of  strains  x  ,  polarization  P,  and  temperature  T 
as  t 


A(x,P,T)  -  \  c"  <*2  +  y2y  +  )  +  ciVyyZz  +  zzxx  +  Vy) 


+  4  cP  (x2  +  yl  +  z*)  +  A  V"  [pz  +  +  Pz]  +  g  rx  pz  +  y  Pz  +  z  Pz] 

2  44 v  x  a  xy  2  A  L  x  y  zJ  S11  x  x  yy  y  z  z 


A  ..H 


+  g10[x  (P2  +  P2)  +  z  (P2  +  P2)  -I-  y  (P2  +  P2)] 
612L  xv  x  z'  z  v  x  y'  Jy  z  x/J 


+  g..[y  P  P  +zPP  +  x  P  P  ], 
“44  J  z  yz  xzx  yyxJ’ 


(1) 


plus  terms  of  higher  order  in  P.  For  a  ferroelectric  material,  x"  is 
a  linear  function  of  temperature  T;  c^,  etc,  are  elastic  stiffness 
constants  at  constant  polarization;  Pz,  etc,  are  components  of  the 
polarization;  zl,  etc,  are  strains;  and  gi^,  etc,  are  electrostrictive 
constants.  The  BaTiO^  ceramic  considered  here  is  macroscopically 
isotropic  material  strained  by  the  shock  pressure  in  the  z-direction. 
The  free  energy  is  a  special  case  of  the  free  energy  (eq  1).  It  is 

A(*,P,T)  .  A  X"  P l  +  A  +  g^P*.  (2) 

We  refer  the  free  energy  to  the  zero  strain,  zero  polarization  state. 
The  differential  free  energy  is 


dA  =  -SdT  -  Z  da  +  E  dP  ,  (3) 

z  z  z  z 

S  is  the  entropy,  Zz  a  stress  component,  and  Ez  an  electric-field 
component. 
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as 


Now  the  inverse  susceptibility  at  constant  strain  is  defined 


z  _ 
X  = 


z 

z 


then  gives  for  constant 


(4) 


the  same.  In  cgs  units, 

z  4tt  4it 

x  »  ~  ~ • 

e  is  the  relative  dielectric  constant,  c  -  1  is  approximately  equal 
to  e,  since  for  BaTiO_, e>>l .  The  slope  of  xZ  with  respect  to  z?  is 
twice  the  9 


Use  of  the  Legendre  transformation  (eq  3) 
temperature  stress 


z 


2 

and  we  see  from  (eq  2)  that  x  and  x"  are 


dxl  _  d _  /92A\ 

dzz  '  dzz  W 


(5) 


A  similar  expression  has  been  obtained  for  the  slope  of  the 
inverse  dielectric  susceptibility  with  respect  to  hydrostatic 
pressure,  p  (ref  9,  16-19).  In  this  case, 

»£  ,  2  Ai  *  *0  ■  (6) 

9p  V 

where 


X 


P  = 


3Bz) 


and  Q..  and  Q._  are  again  electrostrictive  constants.  This  result 
is  derived  most  easily  from  the  Gibbs  elastic  function  G(X,P,T) 

X  represents  stress  components. 

The  g-electrostrictive  constants  relate  the  stress  and  square 
of  the  polarization  of  a  strain-free  material.  The  Q-electrostrictive 
constants  relate  strain  to  the  square  of  the  polarization.  This  can 
be  seen  from  equations  (1)  and  (4).  The  stress  component 


at  zero  strain  is, 


9  A  r  _2  ,  2  _ 2 -t 

9l - [gllPz  +  g12Py  +  Sl2Px- 


(7) 


At  zero  stress  the  strain  component,  z^,  is 


z 

z 


=  QllPz 


+  Q,,P2  + 

12  X 


(8) 


These  are  equations  defining  the  electrostrictive  constants.  The 
Q-constants  may  be  determined  directly  by  measuring  strains  produced 
by  an  induced  polarization.  If, 


then 


P  «,  P 

y  x 


0, 


z 

^1 1=  p2  » 

z 

and,  if  P  =  0,  P  =  P  , 
z  x  y 

then 


z 


x 


(9) 

(10) 


The  Q  constants  are  also  related  to  the  inverse  susceptibility 
versus  pressure  slope  by  expression  (6).  The  Q  and  g  constants  are 
related  by  the  elastic  compliances.  For  a  cubic  crystal  and  also 
for  isotropic  material,  the  relations  derived  from  the  expression 
for  the  free  energy  are  (r.ef  J5,20) 

^11  "  Q12  =  \S11  +  2s12^  (^11  +  2S12) 

and 

(ID 

^11  +  2^1  'i  "  "(Sll  +  2s12)  (811  +  2g12)' 


The  Q  electrostrictive  constants  have  been  measured  directly 
both  for  ceramic  and  single  crystals  (Ref  21,  22)  using  relations 
(9)  and  (10).  The  g^constant  can  now  be  obtained  from  Q  constants 
using  expression  (11).  The  g^  constant  is  related  to  the  slope  of 
inverse  dielectric  constant  versus  uniaxial  strain  by  expression 
(5). 
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In  summary:  Tin*  slope  of  the  Inverse  dielectric  constant  versus 
hvdrost at  ic  pressure  mul  ihe  slope  ol‘  the  same  quantity  versus 
unlaxlnl  strain  can  he  determined  for  a  paraelectrlc  Isotropic 
material  from  the  th.  ami  t).  ,  elect  rust  r  l  et  I  ve  constants  ami  the 
Sjj  ami  the  Sj,  elastic  compTlam.es. 

It  won  hi  also  t  ,<  useful  If  the  actual  transition  strain  anil 
the  actual  transition  shock  pressure  could  also  he  predicted,  for 
a  first  order  transition  the  Inverse  dlectrlc  constant  reaches  a 
minimum  at  the  transition  with  the  value  Itself  discontinuous  there, 
for  a  second  order  transition  the  inverse  dielectric,  constant  Is 
attain  a  minimum  at  the  phase  transit* on,  vanishing  completely  at  that 
point.  There  Is,  however,  no  discontinuity  In  the  value  as  the 
transition  pressure  or  strain  is  crossed,  it  Is  much  easier  to  deal 
with  the  case  of  a  second  order  transition.  Devonshire  and  other* 
have  derived  a  relation  predicting  the  transition  pressure  for  the 
case  of  hydrostatic  pressure,  assuming  that  the  ferroelectric- 
paraelectrlc  transition  could  he  approximated  as  a  second  order  one 
although  the  transition  lr  BaTlQ.  at  atmosphere  pressure  U  a  first 
order  one.  There  Is  fair  agreement  with  experimental  results 
(ref  16- 1 H) .  The  Devonshire  relation  is  derived  most  easily  from 
Gibbs  elastic  function.  A  similar  relation  for  uniaxial  strain 
can  alio  be  derived.  In  this  case  it  Is  most  convenient  to  use  the 
Helmholtz  free  energy.  At  zoro  strain  the  Curle-Welss  law  results 
from  oquation  (2)  if 


(T 


V  • 


(12) 


Tc  is  the  Curie  temperature,  from  equation  (2),  the  inverse 
susceptibility  ut  constant  strain  is 


z 

x  ■ 


*o  (T  -  Tc> . 


(13) 


The  Curie-Meiss  law  is  usually  written  In  terms  of  the  dielectric 
constant  c  nnd  the  Curie  constant  C  as 


L 


(14) 


where 

4n 

c  ■  — — 
z 

X 


12 


mm niiiim  ^ 


hpiat  inn  (i.M  t*  the  Oivlf  Weiss  expression  with  defined  equal  to 
V>  i  ,  It  should  he  pointed  out  that  pquat  Ion  In  valid  only  If 

terms  ot'  h t order  thou  I*'  tit  equal  Ion  tl)  can  ho  Ignored, 

loom  equation  (21,  wo  now  ohtntn  t  ho  expression 


1 

V 

i 

m 

-c**>  •% 

r  "  V  '  J<n*, 

tm 

* 

>\ « 

for  a  second  order  trails  i t ion, 
nt  the  t  rails  it  to  , 

t  •  h  at  the  ?  rails  it  ion 

1 heref ore , 

r  i, 

*  „ 

(16) 

It'  wo  amount)  that  the  shift  ot'  Ourle  temperature  with  strain  U 
Unoat'  over  the  range  between  the  Curie  tempo  tat  tiro  and  t'oom 
temperature,  equation  (16)  becomes  a  relation,  giving  the  shift 
of  Curie  temperature  with  strain. 

1 .  11IKUHY  }W  MKA8URKHKNT 

The  dielectric  constant  of  the  shock-compressed  region  which 
is  in  a  state  of  uniaxial  strain  was  obtained  from  records  of  current 
versus  time.  The  current  is  produced  by  the  depolar laing  effect  of 
the  shock  wave  as  it  propagates  through  the  thickness  of  the  HaTlO] 
disk.  The  current  moves  through  an  effective  short  circuit,  connecting 
electrodes  on  the  disk  surfaces.  The  disk,  showing  electrode  arrange¬ 
ment  and  current  path,  is  shown  in  figure  i.  The  use  of  two  concentric 
electrodes  in  a  guard-ring  arrangement  eliminates  rarefaction  waves 
from  the  lateral  suriace,  The  shock  front  is  plane  and  parallel  to 
the  electrode  in  the  region  between  the  front  inner  and  the  back 
electrode,  The  front  is  a  discrete  step,  less  than  0.2  mm  thick 
in  5-  and  10-mm  samples.  The  pressure  behind  the  shock  front  remains 
constant  during  the  transit  time  of  an  elastfc  (below  the  IIKl.  shock) 
front.  This  was  verified,  using  a  quarts  shock  gauge  of  the  Graham 
type  (ret  23).  The  initial  polarisation  is  axially  directed.  Figure 
2  shows  the  cross  section  in  detail.  We  assume  that  condltons  behind 
the  front  are  constant  during  the  time  the  shock  transits  tire  thickness 
of  the  disk. 


CURRENT 


Figure  1,  ttarlum  Utanate  ceramic  disk  with  alectrodsa.  (The  disks 
were  3. HI  cm  In  diameter  and  either  0.5  or  1  cm  thick). 


Kxpreaaions  have  been  derived  and  reported  previously  for  the 
current  through  a  short  circuit  or  through  a  reaiative  load  (ref  12, 
24-26).  For  a  single  front  and  no  dependence  of  the  dielectric  constant 
on  the  field,  the  current  through  a  path  shunting  the  electrodes  is 


J»  1  +  S 

o  (1  +  StT7  * 


(17) 


where  6Pg  is  the  change  in  polarisation  produced  by  the  shock,  t 
the  transit  time  of  the  shock  wave  between  electrodes,  and  t  is  the 
normalised  time,  which  is  the  actual  time  divided  by  the  transit  time 
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or  (lie  Nltoi'k  I  rool,  A  (a  Dip  p|pi'(ro«|p  mips,  S  ♦  1  U  Dip  ratio  or 
Hip  rliirtl  ( Nlioi'k"i'oin|ii'PMHP(l  titptPi'ial)  to  t  Iip  loti  to!  (imahnckpd) 
d  I  p  Ip<1  r  l  r  I'onatrtiU,  H«v'«o«p  of  l  Iip  high  il  Ip  1p«*  1 r  1 1*  count  ant  of  iIip 
mu !  P  r  i  M  l  *  HOMk'P)ll  III  i  1  I  I  V  rtllll  iltplpit  I:  If  I'OIIMtrtlll  PIP  A|l|troK  llllrt  I  P  I  V 
0Hl'P)ll  lot'  rt  Ul'IOl’  or  All  WltPIl  I'UM  out  IP  rtl'P  UpPtl  i  Wp  iiow 
I'rilvulMlp  Uip  derivallvp,  with  rpMppt'l  to  Dip  ooiipuI  ,  of  Hip  natural 
logai  1 1  htn  of  Dip  i'oitpiiI  . 


«n  A  I* 


Hit  1  +  k’o(l  4'  s) 

l) 


Di(l  ♦  Si)  I,  <1R) 


FRONT 

ELECTRODE 


€f 


*1 


P» 


Pa>6Ps 


BACK 

ELECTRODE 


figure  *- •  Cross  section  of  disk  in  region  of  inner  electrode  where 
plane  geometry  is  valid.  ef  is  the  dielectric  constant 
of  the  shock-compressed  region;  e*  is  the  initial  dielectric 
constant  of  the  material  before  it  is  compressed,  that  is 
the  dielectric  constant  in  front  of  the  shock  front;  6PS  is 
the  polarization  removed  by  the  shock  compression. 


15 


Sliniu  (til  lonitH  Mt‘(>  Independent  of  i  except  the  last, 


<1 

<i7 


iu 


4-  Si) 


-2S__ 
f  t-  s' 


(19) 


The  der Ivatlve  Junl  calculated  can  bn  obtained  from  the  oueil loucope 
record*  of  deiolarlaatlon  current  ver*u*  time.  Once  it  is  known,  S, 
end  consequently  the  ratio  of  initial  to  final  dielectric  constant/ 
in  alao  known.  The  initial  dielectric  constant  Is  a  known  value; 
therefore,  the  dielectric  constant  in  the  shock-compressed  region  is 
determined . 

4 •  EXPERIMENTS  AND  RESULTS 

The  measurement!  of  dielectric  constant  versus  uniaxial  strain 
were  made  on  ceramic  disks  93-wt  percent  BaTit^,  5-wt  percent  CaTi03. 
They  were  3.81  cm  in  diameter  and  either  0.5  or  1  cm  thick.  The 
electrodes  on  the  flat  surfaces  were  of  vapor-dcpositod  gold.  An 
clastic  (pressure  below  the  HEh)  shock  front  was  produced  in  the  disk 
by  impacting  a  stationary  driving  plate  of  either  load  or  aluminum, 
with  a  gas  gun  accelerated  plate  of  plexiglas  aluminum,  or  a 
combination  of  plexiglass  and  aluminum.  The  disk  was  placed  against 
the  stationary  plate.  The  electrodes  were  connected  through  a  low 
resistance  that  was  effectively  a  short  circuit.  The  arrangement  is 
shown  in  figure  3.  Impact  velocity  was  measured,  using  photocells  and 
intercepted  light  beams.  The  disk  was  polarised  to  about  1  percent 
of  the  saturation  polarisation  of  the  ceramic.  This  minimised  the 
magnitude  of  the  internal  fields  during  the  time  the  current  pulse  was 
produced.  The  internal  fields  that  were  estimated  were  low  enough 
to  avoid  shifts  in  the  dielectric  constant  due  to  field  dependence. 

The  results  thus  give  essentially  the  zero-field  dielectric  constant. 

The  current  pulses  were  displayed  on  an  oscilloscope  and  photographed; 
photographic  records  of  these  current  pulses  are  shown  in  fig  4. 

The  vertical  direction  is  current.  Time  moves  from  left  to  right. 
Superimposed  sine  waves  are  calibration  signals.  The  5.2-kbar  output 
shown  in  figure  4  was  generated  by  the  elastic  precursor  wave  in  a 
1-inch-thick  aluminum  plate  used  In  place  of  the  lead  plate  shown 
in  figure  3.  The  current  pulse  produced  by  the  5.2-kbar  front  is  the 
first  portion  of  the  pulse  in  which  the  slope  is  positive;  the  remaining 
pulse  was  produced  by  a  60-kbar  front. 

The  peculiar  negative  and  positive  leading  spikes  superimposed 
on  the  main  current  pulse,  generally  of  swayback  trapezoidal  form, 
are  due  to  the  depolarization  of  surface  regions  in  the  ceramic  (ref 
14,27-29).  They  are  thought  to  have  u  dielectric  constant  relatively 
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independent  of  the  strain  for  the  values  of  strain  involved  in  the 
measurements  and  contribute  little  to  the  quantity  £  In  i,  except  in 
the  immediate  vicinity  of  the  spikes  themselves.  °T 

Values  of  uniaxial  attain  and  shock  pressure  were  obtained,  using 
standard  procedures  (ref  1)  from  particle  velocity  versus  shock 
velocity  data  found  in  the  literature  (ref  11,  13,  19,  25,  30).  When 
guarded  specimens  were  used,  iL  In  i  was  obtained  from  a  digital 
reduction  of  waveforms  for  ab^St  10  values  of  x  between  0.2  and  0.8. 
These  were  then  averaged  and  the  standard  deviation  of  the  values 
averaged  was  calculated.  There  is  always  some  slight  angle  (about 
4'  of  arc)  between  the  impacting  plates.  This  results  in  a  small  tilt 
on  the  moving  front.  The  effect  of  tilt  on  the  output  current  wave¬ 
form  is  to  produce  a  distortion  of  the  curve  form.  The  analysis  by 
Wittikindt  (ref  24)  shows  that  for  intermediate  values  of  t,  the 
difference  in  the  slope  of  distorted  and  undistorted  current  wave¬ 
forms  was  mininuzed.  For  unguarded  arrangement,  however,  the  results 
are  additionally  distorted  for  the  larger  values  of  x  by  rarefaction 
waves  originating  at  the  lateral  surfaces.  All  unguarded  specimens 
were  0.5  cm  thick;  when  x  =>  0.2,  lateral  rarefactions  were  disturbing 
about  7  percent  of  the  area.  For  this  reason,  only  qualitative  results 
could  be  obtained  from  the  unguarded  specimens.  An  exception  was  the 
unguarded  26-kbar  shot  where  a  relatively  hig^impact  velocity  produced 
a  low  apparent  tilt,  allowing  a  good  value  of  -In  i  to  be  obtained 
from  a  small  value  of  x.  The  results  and  some  experimental  details 
are  summarized  in  Tlblt  1  • 

5.  COMPARISON  OF  THEORY  AND  EXPERIMENT 


There  are  no  direct  measurements  of  g,  j  for  BaTi03  ceramic. 

There  is,  however,  a  direct  measurement  of  Qj  j  due  to  Schmidt  (ref  21) 
who  directly  measured  displacement  versus  polarization  for  a  ceramic 
both  below  and  above  the  Curie  point.  He  obtained  a  value  of 
Qj  1  =  0.6  x  10~  cm  /dyne  below  the  Curie  point  and  0.7  x  10"12 
cm  ^dyne  above  the  Curie  point.  He  observed  a  decrease  in  value  in 
the  immediate  vicinity  of  the  Curie  point.  Estimates  of  the  electro- 
strictive  constants  for  the  ceramic  are  also  available, using  single 
crystal  constants  and„  averaging  formulas  (ref  22).  These  are 
Qll  =  0.7  x  10-12  cm  /'dyne  and  Q[2  “  -0.2  cm2  /dyne.  There  is  also 
an  unpublished  estimate  of  Qj  1  and  Q2  (due  to  Saxe  and  quoted  by 
Schj^-dt ,  gef  21):  Qj  j  =  0.64  x  lO”1^  cm  ?/dyne  and  Q|  2  =  -0.138  x 
10  cm  /dyne. 

We  can  also  obtain  a  value  of  Q12  if  we  know  the  value  of  Q 

and  also  the  slope  of  the  line  of  inverse  susceptibility  versus 
hydrostatic  pressure  in  the  paraelectric  region.  The  value  of  3Xp/3p 
obtained  from  experimental  measurements  of  dielectric  constant 
versus  hydrostatic  pressure  by  Sahara  (ref  10)  is 
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2(Qn  +  2Q12) 


0.45  x  10"12 


cm2 /dyne  • 


(20) 


This  value  used  in  conjunction  with  Schmidt's  value  for  Q,  ^  yields 
Q12>  From  the  literature  (ref  31),  we  obtain  values  of  the  elastic 
compliances 


s1 1  =  8.7  :c  10  13  (cgs) 
Sj2  =  -2.9  x  10  13  (cgs) 


(21) 


From  the  relations  due  to  Devonshire  and  quoted  in  section  2,  we  have 


g 

11 


2  Qi ]  ‘  Q12 

3  sn  -  s12 


+  2QI2 
+  2s12 


Table  II  summarizes  the  various  values  of  Qi j  and  Q12  and  also 
calculated  values  of  gj  j  . 


Table  II.  Values  of  Q  and  Q  2  for  ceramic  barium  titanite 
and  value  of  g.. ..  calculated  from  these  (g„  is  dimensionless 
in  c0s  units)  f1  U 


Q-Ll(cm2/dyne) 

Q12(cm2/dyne) 

Source 

811 

0.7  x  10“12 

-0.2  x  10'12 

| 

Single  crystal  average 
(Jona  and  Shirane) 

-0.86 

0.64  x  10'12 

-17 

-0.14  x  10  11 

Saxe 

-0.86 

0.70  x  10"12 

-0.24  x  10"12 

Schmidt  (above  120°C)  § 
hydrostatic  data 

-0.80 
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The  value  of  gj  1  is  thus  about  -0.84.  3'he  predicted  value  of 


3  2  ^ 

- — —  ■  - — -U-  is - -0.1  3  (cgs)  ,  (22) 

dz  4tt  4  it 


The  value  obtained  by  fitting  a  straight  line  to  the  quantitative 
results  (between  8.5  and  26  kbars)  in  Table  I  is 


"  -0.11  ±  0.04  (cgs)  , 


(23) 


The  experimental  points  and  fitted  line  are  shown  in  figure  5.  The 
results  from  shot  1  and  shot  2  were  considered  inadequate  for  reduc¬ 
tion  to  numerical  values  of  inverse  dielectric  constants  and  are  dis¬ 
cussed  In  the  next  section. 

From  (.eq  16),  the  transition  strain  at  room  temperature  (23°C) 
i s  given  by 


-<T  -  Tc)4tt 
zz  2g1'1C 


(24) 


T^  is  112.5  and  C  ■  1.5  x  10s  (ref  6).  On  the  basis  of  table  II,  we 
estimate  gjj  -  -0.84.  The  transition  strain  at  room  temperature  is  then 


zz  =  -0.0044  , 


The  elastic  shock  velocity  in  BaTi03  is  approximately  6  x  105 
cm/sec  and  the  initial  density  (p0)  is  5.5  gm/cm3.  Using  Hugoniot 
relations 


^shock 


0o 


u 

s  p 


(25) 


where  u8  and  Up  are  shock  and  particle  velocity  and 


z 

z 


u 

=*  -  _E 
us’ 


(26) 


we  obtain  the  results  that  at  the  transition  the  particle  velocity 
Up  is  2600  cm/sec  and  the  shock  pressure,  8.6  kbarn. 
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6. 


CONCLUSIONS  AND  DISCUSSION 


The  agreement  between  the  experimental  and  calculated  value  of 
the  slope  of  inverse  dielectric  constant  versus  strain  for  the  para- 
electric  region  suggests  that  in  the  measured  region  between  8.5  and 
26  kbars  the  ceramic  is  in  a  paraelectric  state, and  that  the  transition, 
therefore,  occurs  at  a  shock  pressure  equal  to  or  below  8.5  kbars. 

This  is  consistent  with  the  calculated  shock  transition  pressure  of  about 
8  kbar  and  also  the  value  inferred  by  Doran  (ref  11)  on  the  basis  of  his 
experiments  and  also  those  of  Binder.  The  consistency  between  the  theory  and 
what  are  now  considerable  experimental  results  support  the  contention 
that  shock-induced  ferroelectric-paraelectric  transition  occurs  at 
about  8  kbars  at  room  temperature.  This  value  is  quite  different  from 
the  value  of  hydrostatic  pressure  producing  a  transition-  A  relation 
predicting  the  hydrostatic  transition  pressure  at  room  temperature 
(ref  16-19)  is 


CT  -  T„)4tt 


2(Qn  +  2Q12)c 

It  is  analogous  to  the  relation  predicting  the  transition  uniaxial 
strain.  The  relation  yields  the  results  p  -  17  kbars  where  the  value 
of  2(Q11  +  2Q12)  is  that  obtained  from  the  slope  8-1/3  p  for  pressures 
just  above  the  transition  pressure.  The  slope  was  obtained  from  experi¬ 
mental  measurements  of  Samara  (ref  10) .  This  result  is  in  fairly  good 
agreement  with  the  transition  pressure  20.1  +  0.2  (ref  9)  obtained 
from  measurements  of  four  samples  of  ceramic  BaTi03  (5-wt  percent 
CaTiOg)  and  three  samples  of  technical  grade  BaTi03,  ceramic.  It 
should  be  emphasized  that  the  effect  of  the  calcium  additive  on  the 
transition  pressure  was  small — the  two  materials  exhibiting  essentially 
the  same  20-kbar  transition  pressure. 

It  should  be  mentioned  here  that  the  linear  relation  used  to 
calculate  the  strain  has  been  assumed  to  hold  over  the  full  range 
between  the  room  temperature  and  the  transition  temperature.  This 
assumption  has  been  experimentally  verified  for  the  case  of  hydrostatic 
pressure  (ref  6) . 

Examination  of  t  as  a  function  of  hydrostatic  pressure  (ref  6-9) 
indicates  clearly  that  at  the  transition  pressure  the  value  of  1/e 
is  relatively  large.  The  results  of  the  shock-compression  measurements 
as  seen  in  figure  5  appear  inconsistent  with  the  assumption  l/e-+0  at 
the  transition.  This  inconsistency  can  be  explained  by  assuming 
that  the  measured  capacitance  at  a  transition  is  the  value  of  non- 
ferroelectric  capacitance  layer  electrically  in  series  with  the  capacity 
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"I  t  Itp  Imi  1 !•  i'I  till'  mill  pi  lull  ll  I*  ill*"  I'lmnll'ln  tlirtl  I  Iip  Ihvpi  In 
ft  •h>Ii*iIi|i  Imi  w  1 1  It  I  *  mm  1 1  I  ini  <  Imi  Mi  I  pi  1  nl  1 1  illftrtPiil  limn  llwil 
i'I  I  Iip  I'lill.,  Hip  >  n|nii  llv  •*»  **ti<  It  lavPlw  wmilil  I'M  ipIuIIvpIv 
I  IIiIp|ipii>IpiiI  iif  |i|Pn*uip,  limn,  I  Iip  nipmmipil  value  nf  *'•  ••  /•!# 

wmilil  Ho  I  •  limine.  Hip  i  til  I  pi  I  VrtliiP  fm  tin*  I  I  Mint  1 1  I  mi  yiesmiie 
ipiilil  I  lii'ii  Iip  ini  I  pi  llv  i  rt  1 1  ii  I  «l  pil  |*y  .1  nr  urn  I  ii||  I/'  »  IK  lV*m|hle 
source*  nf  nihil  Inveis  are  nuf  hard  In  fliiil,  Surface  lavris  with 
|irn|ipi  l  I pn  ill  ffei  piiI  f  rom  llml  of  flip  ImlK  have  Imon  ii|u«p|\pd  Imtli 
In  ceramic  mill  single  crystal*  (ref  ,t;i 

I  Imi  f||P  i’pn  I  ilmi  I  Vrtlllr  of  I  ho  Inverse  dlclcrtl'ii  constant  at 
flip  transition  results  fi'oin  ah  extraneous  effect  *  further  supported 
hy  flip  fact  flirt!  flip  VrtliiP  of  l/«  «f  the  transition  when  produced  l»v 
hvdrost  «f  i 0  pressure,  is  extremely  sensitive  to  il»p  particular  type 
of  sample ,  Very  different  value*  of  l/»  are  obtained  for  ceramic  with 
!*  i*t  pp  repot  I'rtilO^  rtilil  1 1  I  vp  ,  on  up  nt  i  rt  t  ly  |UH'p  ceramic  mill  Mingle 
crystals  (t'pf  i>) ,  The  value*  of  flip  slope*  rtilil  transition  pressure*  are, 
however,  approximately  flip  same  am  would  lip  expected  in  flip  presence 
of  h  surface  layer  whose  dielectric  const  ant  U  Independent  of  pressure, 

It  would  littvo  hpoii  vrthirthlp  to  determine  thp  actual  minimum  in 
thp  curve  of  invprMp  dielectric  constant  versus  strain  hy  making 
quant  itrtt  Ivp  measurement*  In  tho  ranye  between  0  and  M  Khar*  ■  Iip  low  thp 
imp  I  toil  U'wiimI  i  ton,  Measurement*  in  thin  you  ion  are  difficult,  because 
low- impact  velocities  ttf'P  required  to  produce  low  strains  rtiul  iow- 
impact  vp  I  oc  i  t  i  pm  produce  fronts  with  htjjh  tilts  from  which  it  is 
il iff! cult  to  obtain  yooil  value*  of  Inverse  dielectric  constant.  The 
result#  ohtalnml  in  shot  l  aiul  shot  2  were  thought  to  1>p  too  ills- 
tortpil  rpltttive  to  pulsp  lenyth  far  reduction  to  quantitative 
values  of  inverse  dielectric  constant.  The  positive  slopes  of  the 
wave  form  however  clearly  Indicate  that  the  inverse  dielectric 
constant  at  these  points  is  greater  than  the  room  temperature  atmos¬ 
pheric  pressure  value  of  8,5  x  10  .  This  is  consistent  with  a 

minimum  in  the  curve  in  the  vicinity  of  8  Khars, 
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Measurements  on  BaTiO..  ceramic  were  made  of  inverse  dielectric 
constant  versus  uniaxial  strain  for  values  of  strain  ranging  between 
h.0043  and  -0,013.  The  strain  was  produced  in  the  region  behind  a 
shock  front  generated  in  a  slab  of  the  material  by  impacting  it  with  a 
flying  plate.  The  plate  was  accelerated  in  a  gas  gun.  The  slabs  were 
disks  with  electrodes  on  the  flat  surfaces.  The  normal  stress  produc¬ 
ing  these  strains  ranged  between  8.5  and  26  kbars.  The  dielectric  con¬ 
stant  was  obtained  from  oscillograms  of  the  shock-induced  depolariza¬ 
tion  current  through  a  circuit  connecting  the  electrodes.  A  straight 
line  was  fitted  to  the  curve  of  inverse  dielectric  constant  versus 
strain.  The  slope  of  the  line  was -0.11  ±  0.04.  Thermodynamic  theory 
was  used  to  calculate  this  slope  assuming  the  material  was  in  the  para- 
electric  phase  for  the  range  of  experimental  shock  pressures.  The 
calculated  value  was -0.13.  The  calculated  value  is  in  agreement  with 
the  experimental  value  and  the  assumption  that  the  material  is  in  the 
paraelectric  phase  for  shock  pressures  between  8.5  and  26  kbars.  The 
agreement  implies  that  the  transition  at  approximately  room  temperature 
is  produced  by  a  shock  pressure  less  than  or  equal  to  about  8.5  kbars - 
This  is  markedly  different  from  the  experimental  lesults  for  hydro¬ 
static  pressure  where  the  transition  occurs  at  20  kbars  at  room 
temperature,  and  the  material  is  in  its  paraelectric  phase  above  20  kbars. 
In  both  cases  the  temperature  is  approximately  room  temperature.  The 
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transition  strain  and  shock  pressure  was  calculated  using  thermodynamic 
theory  and  assuming  that  the  transition  could  be  approximated  as  a 
second  order  one.  The  result  was  a  transition  strain  of  -0.0044  and  a 
transition  shock  pressure  of  8.6  kbars.  This  is  consistent  with  the 
preceding  results  and  also  supported  by  experimental  findings  by  Doran 
who  observed  a  discontinuity  in  sound  velocity  at  that  shock  pressure. 
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